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ABSTRACT 

We present Hubble Space Telescope optical and ultraviolet photometry for five nearby, thermally 
emitting neutron stars. With these measurements, all seven such objects have confirmed optical 
and ultraviolet counterparts. Combining our data with archival space-based photometry, we present 
spectral energy distributions for all sources and measure the "optical excess" : the factor by which the 
measured photometry exceeds that extrapolated from X-ray spectra. We find that the majority have 
optical and ultraviolet fluxes that are inconsistent with that expected from thermal (Rayleigh- Jeans) 
emission, exhibiting more flux at longer wavelengths. We also find that most objects have optical 
excesses between 5 and 12, but that one object (RX J2143. 0+0654) exceeds the X-ray extrapolation 
by a factor of more than 50 at 5000 A, and that this is robust to uncertainties in the X-ray spectra 
and absorption. We consider explanations for this ranging from atmospheric effects, magnetospheric 
emission, and resonant scattering, but find that none is satisfactory. 

Subject headings: Stars: Pulsars: Individual: Alphanumeric: (RX J0420.0-5022, RX J0806.4-4I23, 
RX J0720.4-3125, RX J1308.6+2127, RX JI605.3+3249, RX J1856. 5-3754, 
RX J2I43.0+0654)— Stars: Neutron— X-Rays: Stars 



1. INTRODUCTION 

Among the nearby neutron stars known, seven show 
emission that appears predominantly thermal Q with in- 
ferred temperatures of ~ I0 6 K. These so-called isolated 
neutron stars (INSs) are young, < I Myr old, and the 
thermal emission is thought to be due to residual heat 
as their X-ray luminosities are considerably more than 
their spin-down luminosities E. They differ from simi- 
larly aged pulsars not only in the absence of non-thermal 
radio and X-ray emission, but also in their long, 3- 
10 s spin pe riods and la r ge, ~ 10 13 G magnetic fields 
(for reviews, IHaberl|[2007l:lvan Kerkwiik fc Kapladl20"07t 
lKaplanll2008l) . 

The INSs have attracted much attention, in part be- 
cause of the hope that their properties could be used to 
constrain the poorly understood behavior of matter in 
their ultra-dense interiors: since th e details of the neu- 
tron star's interior affect its radius (|Lattimer fc Pra kash 
I2007H and cooling history (|Yakovlev fc Pethickl I2004D . 
the wide range of theoretical possibilities can be lim- 
ited with data. In this respect, the thermal emission 
from INSs is particularly interesting, since one might de- 
rive constraints on mass, radius, and cooling history from 
spectral properties such as effective temperature, surface 
gravity, and gravitational redshift. 

Progress has been stymied, however, by difficulties in 
interpreting the observed spectra: at present, both the 
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composition and state of matter in the photosphere re- 
main unknown, with hydrogen, helium, and mid-Z ele- 
ments (O, N, Ne, etc.) in states ranging from gaseous to 
condensed all being consider ed (see, e.g., contributions to 
iPage. Turolla. fc Zandl2006l) . Below, we use the bright- 
est and best studied INS to illustrate those problems and 
the part played by optical and ultraviolet measurements. 

1.1. The Puzzle of RX J 1856. 5- 3754 

The emission of INSs was known to be roughly ther- 
mal, but it came as a major surprise when long Chandra 
and XMM observations found that the spectrum of the 
brightest INS, RX J1856. 5-3754, could be reproduced 
with a featureless black body (e.g.. IBurwitz et al~ll2001l 
2003). The reason for the surprise was that, for a light 
element (H or He) atmosphere, the spectrum may be fea- 
tureless but the high-energy tail should be harder than 
a Wien tail (because opacity decreases towards higher 
energies and hence one should see deeper, hotter lay- 
ers), while for heavier elements, the overall shape may be 
blackbody-like but spectral features should be present. 
When combined with the p arallax deter- 

(|Walterl [2001 



mined from HST observations 
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Walter fc Lattimerl 1200% I Walter et all 120101) , the 
implied blackbody radius was too small for a neutron 
star, leading to speculation that the object might be 
a "quark star" (|Drake et al.l 120021 ). The suggestion of 
a quark star, however, ignored that the optical and 
ultraviolet measurements were inconsistent with a single 
blackbody, being in excess by a factor ~ 8 over the 
extrapolation of the X-rays (IWalter fc Matthews! Il997t 
Ivan Kerkwiik fc Kulkarnl I2001L and see Fig. [Tj). This 
excess has been taken as evidence for two regions on 
the surface: a hotter one primarily responsible for the 
X-ray emission, and a larger, cooler on e responsible for 
the optical (e.g. JBraie fc Ro mani 2002). In the context 
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Table 1 

Summary of HST optical and UV observations 



RX J Instrument and Filter Date Exp. ST Mag. a 

CUT) (s) 



For the hydrogen models described above - and gener- 
ally for atmospheric models - the optical excess depends 
on the magnetic field strength. For RX J1856.5— 3754, 
the field could only be constrained by the absence of 
X-ray absorption features, but for other sources the sit- 
uation is better: these have absorption features, which 
any model will have to reproduce simultaneously with 
the optical excess, and at field strengths consistent with 
those derived from timing. 

For all INSs, optical counterparts have been searched 
for in deep optical observations, leading to four se- 
cure identifications and two likely identifications. So 
far, apart from RX J1856.5— 3754, only for the second- 
brightest INS, RX J0720.4-3125, has the full optical 
to UV range been probed. The interpretation is dif- 
ficult, though, partly because the fluxes are somewhat 
inconsistent with tho s e of a Rayleigh- Je ans spectrum 
([Kaplan et all I2003bb iMotch et all 120031 ) . and partly 
because, uniquely among the INSs, its X-ray spec- 
trum changed (which poses problems and possibilities 
all of its own: Ide Vries et all [2001 lHaberl et all [200l 
Ivan Kerkwijk et al.ll2007HHohle et al.ll2009l ). 

The possible presence of non-thermal emission in 
the optical, also suggested f or R X J1605. 3+3249 
(|Motch et all 120051: IZane et al.l I2006D . makes it diffi- 
cult to interpret the optical emission. However, for 
RX J0720.4-3125 and RX J1856. 5-3754 the ultravio- 
let emission does seem completely thermal. This is sim- 
ilar to what is found fo r comparably-aged radio pulsars 
such as PSR B0 656+14 f Pavlov. Weltv. k Cordovall997l: 
Shibanov et al.l 120051 ) and Geminga ([Kargaltsev et al.l 
20051) . (Note, though, that for those pulsars the ultravi- 
olet emission is generally consistent with the X-ray ex- 
trapolation.) 

Given the above, we obtained HST optical and UV 
photometry for those five INSs that had not been studied 
in detail before. We complemented this with other HST 
photometry for the INSs (reanalyzing the data where 
necessary); we restrict our analysis to HST data both 
for the high signal-to-noise that it often implies as well 
as the uniform quality of the calibration. 

In Section [2] we present our new data, and give de- 
tails concerning the identification and photometry of the 
counterparts. In Section [31 we compare our results to 
previous (usually ground-based) searches, and perform 
detailed spectral fitting of the optical/UV data for all 7 
INSs, with and without reference to the X-ray spectra. 
We present our discussion and conclusions in Section |4l 

2. OBSERVATIONS AND ANALYSIS 

Of the seven confirmed INSs discovered by ROSAT, 
three do not have confirmed optical counterparts: for 
RX J0420.0-5022 and RX J1308.6+2127, the associa- 
tions are based only on positional coincidence, without 
confirmation based on spectrum or proper motion. Fur- 
thermore, RX J1605. 3+3249 and RX J2143.0+0654 did 
not have ultraviolet photometry. We observed all of these 
with HST in the optical and ultraviolet. A log of the ob- 
servations and our photometry are given in Table [TJ 

2.1. Image Reduction and Combination 

The data for each of the five targets consist of mul- 
tiple exposures covering the SDSS g' band (F475W fil- 
ter, centered near 4700 A) with the Advanced Camera 



0420.0-5022 ACS/WFC F475W 
ACS/SBC F140LP 

0806.4-4123 ACS/WFC F475W 
ACS/SBC F140LP 

1308.6+2127 ACS/WFC F475W 
ACS/SBC F140LP 

1605.3+3249 ACS/WFC F606W 
ACS/WFC F475W 
ACS/SBC F140LP 

2143.0+0654 ACS/WFC F475W 
ACS/SBC F140LP 



2009-10-24 2424 27.85 + 0.25 

2009- 12-16 2856 22.85 + 0.08 

2010- 05-19 4868 27.92 + 0.22 
2009-12-15 5692 23.61 + 0.11 
2009-08-08 4676 27.97 + 0.25 
2009-08-06 5504 23.79 + 0.09 
2005-02-06 4728 27.26 + 0.07 
2009-08-01 2256 26.61 ± 0.09 

2009- 01-18 2688 22.62 + 0.07 

2010- 05-19 7076 26.31 ± 0.05 
2010-05-23 8296 23.38 + 0.06 



a Values are magnitudes in the STMAG system that have been 
corrected for finite apertures. Aperture corrections amounted 
to 0.18mag for the F475W/F606W data and 0.33mag for the 
F140LP data. We estimate an additional 5% systematic un- 
certainty on the SBC/F140LP photometry owing to uncertain- 
ties in the aperture correction. The zero points used were 
26.67, 25.75, and 20.316 for F606W, F475W, and F140LP, re- 
spectively, taken from the revised calibration for AC S given at 
http : //www. stsci . edu/hst/acs/analysis/zeropoints/ 

of this model, the radius is not small but rather puz- 
zlingly large. Furthermore, the a mplitude of the X-ray 
pulsa tions is surprisingly tiny (jTiengo k M ercghetti 

As an alternative, it has been suggested that 
the surface is condensed, but covered by a thin 
gaseous layer of hydrogen (s e e IHo et all 120071; also 
Motch. Zavlin. k Haberlll2003l: IZane. Turolla. k Drake! 
20041 ). Hydrogen can unify the optical and X-ray emis- 
sion, since compared to a blackbody, a hydrogen at- 
mosphere will appear to show an optical excess by an 
amount that depends on the m agnetic field strength 
(|Ho. Potekhin. k Chabrierl200 8[) , and the pulsations will 
be determ ined by non-uniformities in the temperature 
pol 120071 ). For a thick hydrogen atmosphere one ex- 
pects a high-energy tail, but a suitably thin atmosphere 
will be transparent at high energies, allowing one to see 
the blackbody-like emission from the condensed surface 
below. 

By construction, the above model reproduces the 
spectrum. But it has other advantages: (i) It re- 
solves the problem of the small radius: because of 
its non-gray opacities, the temperature is smaller and 
the radius larger than for a pure blackbody model; 
ii) For the dipolar magnet ic field of 1.5 x 10 13 G 
van Kerkwiik k Kaplanll2008f) . a heavy-elemen t surface 
could indeed be condensed (jMedin k Lail 120071 ) . though 
we note that the field is stronger than inferred by 
IHo et al.l ([20071 ): (iii) For other INSs, hydrogen might be 
responsible for the observed spectral features (see lHaberU 
[20071 Ivan Kerkwiik k Kaplan! [20071 ): and (iv) The ap- 
pearance of a thin layer of hydrogen may be the easi- 
est explanation for the change in spectrum o bserved for 
RX J0720.4-3125 (jvan Kerkwiik et al.ll2007l ). 

While promising, this model arguably is contrived, and 
certainly for RX J1856.5— 3754 it leaves it unclear how to 
proceed to test it. Fortunately, however, the other INSs 
provide much more information. 

1.2. Expectations for the other Isolated Neutron Stars 
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Figure 1. X-ray to optical SEDs of RX J1856.5-3754 (left) and RX J2143.0+0654 (right). We plot archival X-ray data fro m XMM (blue) 
and CXO (green) along with the best-fit models (dashed red lines); the model for RX J1856.5— 3754 is just a blackbody (Burw itz et al.l 
2003) while that for RX J2143. 0+0654 includes a broad absorption line (Kaplan & van Kcrkwijk 2009b). The unabsorbed models are the 
solid blue lines. In the optical/UV we plot the HST data from van Kcrkwiik & Kulkarni (2001) and this paper along with the best-fit 
power-law: the absorbed power-law is the black dashed line, while the unabsorbed power-law is the solid black line. 



for Surveys (ACS) Wide Field Channel (WFC) and the 
near-UV (F140LP filter, centered near f500A) with the 
ACS Solar Blind Channel (SBC). For both the WFC and 
SBC data we used four sub-pixel dithered exposures per 
orbit and 1-3 orbits per object; for the WFC the object 
was placed at the center of the WFC1 detector so that we 
did not have to deal with the gap between the detectors. 
For completeness, we also re-analyzed the ACS/WFC 
F606W data on RX J1605. 3+3249, si nce no photometr y 
is given in the original publication bv lZane 

WFC images obtained after Servicing Mission 4 (SM4) 
show charge trails along the CCD columns and faint 
stripes along the rows. The charge trails are due to 
the CCDs' degrading charge transfer efficiency (CTE), 
while the stripes reflect a problem in the new electron- 
ics. These artifacts can significantly affect the detection 
of very faint sources as well as photometry and astrom- 
etry. We used the tasks acs_destripe (version 0.2.1) 
and CteCorr (also version 0.2.1, with the reference file 
pctef ile_101109.f its; see lAnderson fc Bedm1l2010D to 
correct the images, finding noticeably improvements in 
appearance. Checking the resulting photometry against 
that of raw images, we confirmed the conclusions of 
lAnderson fe Bedinl (|20100 that their algorithm restores 
photometric accuracy. Compa ring against the averag e 
CTE corrections computed by IChiaberge et al.l (|2009f ). 
we found good agreement, although with some scatter 
(agai n, similar to what was found by Anders on fc Bedinl 
120101 ). For the F606W data on RX J1605.3+3249, taken 
before SM4, the stripes are not present and we could not 
use the pixel-based CTE correction since the data were 
taken before SM4 and would require different calibration 
parameters, but other methods are possible to correct 
the reduced CTE losses. The SBC data do not suffer 
from such artifacts. 

For each source in Table [TJ 'destriped' and CTE- 
corrected (as necessary) images wer e drizzled using 
multidrizzle (jKoekemoer e t al. 2003) onto a single im- 
age. We experimented with the drizzle parameter^ to 

5 Following http: / /www. stsci . edu/hst/HST_overview/document 



balance resolution and uniformity of sampling, settling 
on a set that gave good image quality: for the WFC 
data we used a final pixel scale of 0.04"pixel _1 and a 
pixel fraction of 0.9, while for SBC we used a pixel scale 
of Q.03"pixel _ and a pixel fraction of 0.9. We verified 
that our photometry did not depend on the choices that 
we made. 

2.2. Astrometry and Counterpart Identification 

For all INSs, the SBC images show an obvious, rather 
bright counterpart, with at best a few other faint sources 
present. In contrast, in the WFC images, the coun- 
terparts are faint and there are many other, brighter 
stars. Previous identification of optical counterparts to 
the INSs has relied on absolute astrometry using the 
positions measured in X-rays and approximate ties be- 
tween the X-ray an d optical coordinate systems (e.g., 
iKaplan et afll2003al ). This can lead to uncertainties, as 
the absolute X-ray positions are typically accurate to 
only 0'.'6 or so. However, the SBC images allow us to 
conclusively identify the counterparts with much higher 
accuracy. To do so we started by registering the SBC 
and WFC images. In two cases, we estimated offsets 
from single, bright optical stars that had faint ultravio- 
let counterparts, while in the others we relied on diffuse 
emission from galaxies (smoothing the SBC and WFC 
data as necessary to isolate a single bright component). 

In all cases the registration was unambiguous, and led 
to clear identifications of the optical counterparts to the 
INSs. The SBC and WFC positions matched to within 
< 0'.'04 (since they were taken at most 7 months apart, 
a typical proper motion of 100 mas yr -1 would lead to a 
comparable intrinsic offset). 

For all of the WF C images, we then regi stered the as- 
trometry to 2MASS (|Skrutskie et alJl2006f >. We fit for a 
shift only, leaving the plate-scale and rotation as set by 
the drizzling software. In the case of RX J0806.4— 4123, 
we used 49 stars, restricted to being on the same de- 
tector (WFC1) as the presumed counterpart. For the 
other images, m any fewer stars were available, and we 
s/muife^^g^g^ifra^J^piaill unsaturated and with stellar 
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Figure 2. ACS images of the fields around RX J0420.0-5022 (top left), RX J0806. 4-4123 (top right), RX J1308. 6+2127 (center left), 
RX J1605.3+3249 (center right), and RX J2143.0+0654 (bottom). For each source we plot WFC/F475W and SBC/F140LP images. 
The optical counterparts are labeled with "X"'s; their positio ns agree to better than 0"05 wi th bright sources in the SBC/F140LP 
images. For RX J0420. 0—5022, selected reference sources from Habcrl, Pietsch, & Motch (1999) are labeled. The positions of objects 
"A" and a galaxy to the south-east match between the optical and UV images. For RX J0806.4— 4123, the bright star to the north-east 
(2MASS J08062407— 412223.1) is present on both images with the same position. North is up, and east is to the left, and the images are 
20" X 20". The SBC images have been smoothed with a Gaussian kernel with a = 2 pixels. 



profiles). From the star-to-star variations, we infer that 
our absolute astrometry is accurate to ~ 0'.'2 only, but 
this suffices to show that the positions of the identified 
counterparts agreed with published X-ray positions (and 
generally with previous identifications; see Sect. [3]). We 
defer more accurate relative and absolute astrometry to 
a later paper. 

2.3. Photometry 

We carried out aperture photometry of the final driz- 
zled images using standard procedures in IRAF. For the 
WFC data, we used an aperture of 0'.'2, and measured the 
sky between radii of 0'.'4 and 0'.'8 (for comparison, the ef- 
fective FWHM in ou r drizzled images is O'.'l). Based on 
ISirianni et al.l (|2005D . the aperture correction from 0'.'2 to 
infinite radius for the F475W filter is 0.183 mag. There is 
a small correction to our aperture correction, since some 
of the light from the stellar PSF will appear in the sky 
annulus. This will bias the sky value upward and hence 
make the source appear fainter than it is. The fraction of 
the stellar light in the sky annulus is 0.04 mag, based on 



ISirianni et al.l (|2005[) . Given the ratio between the areas 
of the source aperture and the sky annulus, the implied 
bias to our photometry is 0.003 mag, considerably less 
than our photometric uncertainty. For F606W, we pro- 
ceeded simil arly, except since we could not use the CTE 
correction of lAnderson fc Be din (20l6J) on those data, we 
applied an empirical correction of 0.034 mag to the pho- 
tomet ry, as inferred from the work of IChiaberge et al.l 

pool . 

For the SBC data our procedure was similar. We es- 
timated the sky background by computing a mean of 
the data between radii of l'/2 and 2'.'4; here, we found 
it necessary to turn off the default clipping of outliers, 
as this removed valid data: the sky is so dark that pix- 
els contain only one or two counts, leading to a highly 
non-Gaussian distribution of sky values. Source photom- 
etry was done within a radius of 0'.'24, and we used an 
aperture corr e ction of 0.33 mag, based on the work of 
IProffitt et al.l (|2003T FI. Given the large aperture correc- 

6 IProffitt et"aH |2003) was specifically concerned with the STIS 
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Table 2 

X-ray Blackbodies for the INSs 



RX J 




kT 


R/d 


Ref. 


Notes 




(xlO 20 cm" 2 ) 


(eV) 


(km kpc -1 ) 






0420.0-5022 


2.02 + 0.71 


45.0 + 2.6 


12.2 ± 0.2 


1 


Normalization is from a blackbody fit. 


0720.0-3125 a 


1.04 + 0.02 


88.4 + 0.4 


15.7 + 0.2 


2 


From 2002-Nov, using XMM/EPIC 


0806.4-4123 


1.7 + 0.2 


87.2 + 1.1 


9.6 + 0.6 


3 


Blackbody with 2 Gaussians 


1308.6+2127 


1.8 + 0.2 


102 + 2 


8 + 2 


i 


Normalization recomputed from CTiandra/LETG 


1605.3+3249 


0.98 + 0.19 


92.6 + 0.8 


12.0 + 0.6 


5 


Blackbody with 1 Gaussian in wavelength, XMM/EP1C 


1856.5-3754 


0.95 ± 0.03 


63.5 + 0.2 


36.7 + 0.8 


6 




2143.0+0654 


2.28 + 0.09 


104.0 + 0.4 


6.2 + 0.1 


7 





References. — 1: IHaberl et all J2004D; 2: IHohle et all d2009D; 3: IKaplan fe van Kerkwiikl <2009blk 4: ISchwope et all i2007D ; 5: 
Ivan Kerkwiik et al.| 1 120041) : 6: IBurwitz etaTI d2003D ; 7: IKaplan fc van Kerkwiikl J2009air 

Note. — The normalization and temperature arc as viewed by an observer at infinity. Where multiple spectral fits arc given in a 
paper, we try to be specific about which we have selected. While the best-fit spectra in many cases include absorption components, 
we only list the relevant blackbody parameters here. 

a The X-ray spectrum of RX J0720.4— 3125 evolves ]de Vries et alJl2004D ; while little if any evolution is present in the optical flux 
(Kaplan et al. 2007). wc attempted to find an X-ray spectrum close in time to the HST observations from Kaplan et al.l I l2003bl) which 
were taken between 2001-July and 2002-Feb. 



tion and the possibility of scattered light beyond 1", we 
include a 5% systematic uncertainty on our SBC pho- 
tometry based on the scatter between the PSFs of the 
different objects. The final photometry are given in Ta- 
ble [TJ All photometry is in the STMAG system, which 
is defined so that a source of constant F\ has a constant 
magnitude: STMAG = -2.51og 10 F A - 21.1 with F x in 
erg cm~ 2 s _1 A -1 . 

3. COUNTERPART IDENTIFICATION AND SPECTRAL 
FITTING 

We show the images for all sources in Figure [2] As dis- 
cussed above, in all cases there are faint optical sources at 
the positions of the ultraviolet sources, whose positions 
also agree with the X-ray positions. Thus, our identifi- 
cations are secure. Below, we first compare our results 
with previous work, and then discuss the ultraviolet and 
optical fluxes, both on their own and in relation to X-ray 
measurements . 



RX J0A20. 0-5022— IHaberl et al.l (j2004l) identified a 
possible counterpart, which is roughly consistent in posi- 
tion. However, at B = 26.6±0.3, it is almost a magnitude 
brighter than our counterpart (using a correction from 
STMAG to Vega-based B-band of m F 47 5 w -B = -0.26, 
computed using synphot for a Rayleigh- Jeans spectrum; 
this should be accurate to b etter than !% ) ■ It m ay be 
that the object identified by IHaberl et al.l f|2004T ) was a 
blend of the true counterpart and another object of sim- 
ilar brightness about 0'.'5 to the North (labeled "Y" in 
Figure [2j although it is hard to see in that image) . Our 
photome try is much more cons istent with the object pro- 
posed bv iMignani et al.l (|2009ft . 

FUV detector. However the ACS/SBC is a copy of that detector, 
and the aperture corrections should be similar (accounting for 
slightly different plate-scales). This approach was recommended 
to us by the STScI help desk. Charge scattering in the SBC 
detector creates a halo extending out to roughly 1" that con- 
tains about 20% of the light, according to the ACS Data Handbook 

(http://www.stsci.edu/hst/acs/documents/handbooks/DataHandbookv^/ja^sc^pg,l^i^]^^ f ll^^)ji data" we use below We SUSPect 
accounting tor the large magnitude ot the aperture correction. 7 ,i j. • a , fi ft, 'ft, l 

While our targets do not have enough counts for very accurate tnat tne discrepancies reflect mostly tne more tnan usual 
profiles out to these radii, a comparison of all of th e stars agrees care required to obtain reliable ground-based measure- 
reasonably well with the aperture corrections from IProffitt et"aLl ments for very faint objects (see, e.g., t he discussion for 
COOl) for a wavelength of 1500 A. Rx J0720.4-3125 in IMotch et aUl2003l) . 



RX J '1 308. 6+2127— IKaplan. Kulkarni. fc van Kerkwiikl 
(20023) identified the same counterpart as we did. Us- 
ing STIS with no filter, which gives a broad band cen- 
tered near 4200 A, they measured msoccD = 28.56±0.13. 
Given msoccD — wf475W = 0.36 mag (again, using 
synphot), this is consistent with our photometry. 

RX J1605. 3+3249 — While the counter part is secure 
based on both color s and proper motion (IKaplan et"aH 
I2003al iMotch et all 120051: iZane et al.l |2006j), there~is 
some quest i on re garding the optical photometry, with 
IZane et al.l (j2006) arguing the it cannot be fit with a 
single power law. We do not confirm this: with our re- 
analysis of the F606W data, the space-based photome- 
try is consistent with a power law (see bel o w). F or the 
ground- based observations of IMotch et al.l ([2005), only 
the i?-band point is consistent with what we measure 
(B = 27.22 ± 0.10 and R = 26.9 ± 0.14, and, using 
synphot, mF475W — B = —0.32 mag and mF606W — R = 
0.23). 

RX J2U 3.0+0654 — We ide ntify the same object as tha t 
found by IZane et all (|2008t) and ISchwope et ail (|2009j) . 
These authors each present S-band measurements, of 
B = 27 .4 ± 0.2 and 26.96 ± 0.20, respectively. These 
are only marginally consist ent; our measur e ment is more 
consistent with the latter. IMignani et al.l (|2011| ) set an 
upper limit of R > 27, which is a bit fainter than 
our expected value of R « 26.6 (from synphot we find 
m475w — R = —0.26 mag). T hese discrepancies m ay be 
from calibration uncertainties (jMignani et al.ll201 ll found 
the two -B-band measurements to be consistent given un- 
certainties in atmospheric extinction correction), but at 
least for this source we cannot yet rule out variability or 
spectral curvature. 

Overall, our identifications our consistent with prior 
ones, but some of the fluxes are discrepant. We do 
not undertake here to resolve this, but note that there 
is no evidence for variability or non-power-law optical- 



6 



Kaplan et al. 



3.1. Power-law Fits 

As in Ivan Kerkwiik fc Kulkarnil (pOOl and sim- 
ilar papers, we tried fitting single power-laws 
to just the optical/UV data for all 7 INSs; 
we included the space-based photometr y for 
RX J1856. 5-3754 (Ivan Kerkwiik fc Kulkarril l200lh 
and for RX J0720.4-3125 (|Kaplan et al.1 l2003b[ ). De- 
tailed fitting of HST data for the INSs presents some 
difficulties, since the zero-point fluxes were determined 
for a spectrum that is flat in F\, while the INSs have 
much steeper spectra. Given the often very wide 
bandpasses used, this leads to substantial changes in 
effective wavelength. 

To circumvent these difficulties, we relied on synphot 
to compute the expected STMAG as a function of ex- 
tinction Ay and power-law index a (with F\ oc A~ Q , 
where a Rayleigh- Jeans spectrum would have a = 4) 
on a fine grid, with steps 0.005 mag in Ay and 0.025 in 
a. Then, for each object, we fitted a power law, opti- 
mizing for a and normalization, but keeping the redden- 
ing Ay fixed to the value implied by the column den- 
sity Ah determined from fits to the X-ray spectrum : 
Ay = Ah/1.79 x 10 21 cnT" 2 (IPredehl fe Schmittl IT9951: 
note that this relation varies for different sight lines, but 
the extinction is small enough that the variations do not 
influence our results qualita tively, as we show below). We 
use th e extinction curve of iCardelli. Clayton, fe Mathisl 
(1989), implemented as gal3 in synphot, and we assume 
E(B—V) — Ay/3.1. We note that we cannot fit uniquely 
for Ay based only on our data, as the amount of redden- 
ing is highly degenerate with the power-law index. The 
exact value of TVh (and hence Ay) that we used depends 
somewhat on the calibration of the X-ray data and on the 
assumed spectral model, and could change in the future. 
We discuss this explicitly below. 

The results of our power-law fits are shown in with 
Table [5] and Figure [3] Our uncertainties and \ 2 val- 
ues include a 5% systematic uncertainty for the F140LP 
data, owing to uncertainties in the aperture correction. 
For those objects with only two measurements we are 
left with a fully constrained fit (0 degrees-of- freedom); 
for the others, % 2 is consistent with a good fit. 

The precise results depend slightly on our choice of 
extinction, but the results do not change qualitatively, 
since the extinctions inferred from the X-ray column 
densities are small (Ay < 0.1 mag, and thus A1400 < 
0.3 mag). Indeed, given the quoted uncertainties on Ah, 
the changes are nearly negligible. Somewhat larger sys- 
tematic changes are possible from changes in the red- 
dening law, the ratio of Ah to Ay, or the assumed X- 
ray spectral model. Therefore, in Table [5] we also give 
da I dAy : the change in the power-law slope with change s 
in extinction (see also Ivan Kerkwiik fe Kulkarnil l200l . 
We note, though, that even a 50% change in Ay would 
generally change a by less than 0.1. 

3.2. Comparison with X-ray Spectra 

One of the main purposes of our work is to establish 
reliable estimates of the "optical excess," the factor by 
which the optical/UV flux of the INS exceed the ex- 
trapolation of the X-ray blackbodies. To do this, we 
have collected the best-fit X-ray spectra for all 7 INSs 
in Table [2] In most cases the best-fit spectra are not 



purely blackbodies, but for the purposes of extrapola- 
tion the blackbody component is sufficient. We note that 
the spectrum o f RX J0720.4-3125 changed during 2003 
(jde Vries et all 120041 ) . While there is no evidence that 
the o ptical flux of RX J0720.4— 3125 changed du ring that 
time (jKaplan. van Kerkwiik. fe Andersonll2007l ) , the two 
different X-ray spectra differ by a factor of 1.54 from each 
other, thus introducing a systematic uncertainty into our 
inferred optical/UV excess. It does not, however, affect 
the shape of the excess, as the HST data are well into 
the Rayleigh- Jeans tail and uncertainties in the extinc- 
tion are small. 

For all of the objects we determine the excess rela- 
tive to expected fluxes calculated by passing the extrap- 
olated, reddened X-ray spectrum through the HST filter 
response curves, using synphot. This avoids issues of un- 
certainties in the effective wavelengths or zero-point flux 
that come from using spectra that differ significantly in 
slope from the calibration spectra. We determined for- 
mal uncertainties by repeating the above for values that 
differ by 1-cr in all parameters from the best-fit spectrum, 
taking as l-er uncertainties on the extrapolated magni- 
tudes the maximum and minimum of these variations. 

We list our results in Tables [3] and |U For completeness, 
we include expected values for instrument/filter/object 
combinations beyond those used here, and also give ref- 
erence photometry for an unabsorbed 10 6 K blackbody as 
well as approximate valu es for the wavelength-dependent 
extinction A\/Ay (see Ivan Kerkwijk fc Kulkarnil l200ll 
for an extended discussion). The extinction values actu- 
ally depend slightly on the value of the overall extinction 
Ay, since this changes the shape of the spectrum and 
hence the integral over the filter response. The effect is 
largest for the widest filters, like STIS/50CCD. However, 
for the range of extinction values considered here (Ah = 
(1 - 4) x 10 20 cm" 2 corresponds to Ay w 0.06 - 0.22), 
the changes are less than 1% in A\/Ay. 

Along with the more easily assessed uncertainties on 
kT and (R/d), a possibly significant systematic uncer- 
tainty in our model is variations in extinction. Similar 
to the analysis above, we explore this issue by consid- 
ering a ±50% range on Ay around the nominal value. 
This results in uncertainties on the extrapolated pho- 
tometry ranging from 0.09 mag (for the shortest wave- 
lengths) to 0.03 mag (for the longest), based on a fiducial 
Ah = 10 20 cm~ 2 . For other values of Ah, the uncertain- 
ties are proportionate. Generally, these additional uncer- 
tainties are smaller than the photometric uncertainties, 
and will not affect our results. The measured excesses 
are given in Table [5] and Figure [3] 

4. DISCUSSION & CONCLUSIONS 

4.1. Comparing the INSs to Each Other 

It is clear that there is a wide variation in both the 
normalization and slope of the optical/UV power-laws 
for the INSs. Co mparing to previous work , we co nfirm 
the conclusion of Ivan Kerkwiik fc Kulkarnil ((2001) that 
RX J1856.5— 3754 has a nearly thermal spectrum, close 
to F A oc A~ 4 , while the spectrum of RX J0720.4-3125 
is somewhat fla t ter, co nsiste nt with wha t was found 
bv Kaplan et~aTI (|2003bD and iMotch etaLI (|200l . The 
other objects extend the range, with RX J0420. 0-5022 
also appearing thermal and RX J2143. 0+0654, with a = 
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Figure 3. Optical/UV SEDs for the 7 INSs. For each object (sorted by increasing kT from top to bottom), we plot \ 4 F\ versus wavelength; 
a Rayleigh- Jeans tail would be flat. The best-fit power-laws with ±ltr uncertainties are shown by the cyan lines. The extrapolations of 
the X-ray blackbodies with ±1<t uncertainties are shown by the magenta lines. The fluxes have been corrected for absorption using the 
nominal values of TVjj from Table [2] The red stars show the wavelength of lOfcT. The different filters are labeled at the bottom. 
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Table 3 

Predicted HST UV Photometry from Extrapolated X-ray Blackbodics 



RX J 


Extrapolated STMAG 






9T1S/FUV ST1S/FUV ACS/SBC WFP02 
F25SrF 2 F140LP F170W 


ST1S/NUV 

F25Qtz 



0420.0-5022 24.87 ± 0.07 25.10 ±0.07 25.30 ± 0.07 25.95 ± 0.07 27.03 ±0.07 

0720.0-3125 23.36 ±0.03 23.61 ± 0.03 23.82 ± 0.03 24.48 ± 0.03 25.57 ±0.03 

0806.4- 4123 24.57 ±0.14 24.81 ±0.14 25.01 ±0.14 25.67 ±0.14 26.76 ±0.14 
1308.6+2127 24.79 ±0.55 25.03 ± 0.55 25.24 ± 0.55 25.89 ± 0.55 26.99 ± 0.55 
1605.3+3249 23.88 + 0.11 24.13 ±0.11 24.34 ±0.11 25.00 ±0.11 26.09 ±0.11 

1856.5- 3754 21.88 + 0.05 22.13 + 0.05 22.34 + 0.05 23.00 + 0.05 24.09 + 0.05 
2143.0+0654 25.41 + 0.03 25.64 + 0.03 25.84 + 0.03 26.50 + 0.03 27.59 + 0.03 
10 b K a WM 2O08 2029 2096 22X15 
A cff a (A) 1355 1442 1517 1769 2283 
A x /A v a 3.11 2.84 2.71 2.65 2.58 



Note. — Wc give the instrument name and filter for all UV observations that we consider; 
not all objects were observed with all combinations, and see Tablc^for the optical observations. 
The 1-a uncertainties arc based on the uncertainties in the X-ray spectra (Tablc[2j and do not 
include any other systematic terms. 

a The STMAGs for a 10 6 K (86 eV) blackbody with A v = and normalized to m F4B ow = 25. 
The wavelengths given are approximate effective wavelengths for the instrument/filter and a 
typical INS spectrum. Wavelength-dependent extinction A\jAy for a 10 6 K (S6eV) black- 
body, calculated for Ay = 0.05. Effective wavelengths and extinctions follow the definition of 
van Kerkwiik fc K ulkarni (2 00iT) . 



Table 4 

Predicted HST Optical Photometry from Extrapolated X-ray Blackbodics 



RX J 






Extrapolated STMAG 








WFP02 
F300W 


WFPC2 

F450W 


ACS/WFC ST1S/OOD ACS/WFO 
F475W 50CCD F606W 


WFPC2 
F606W 


ST1S/CCD 
F28x50LP 



0420.0-5022 
0720.0-3125 

0806.4- 4123 
1308.6+2127 
1605.3+3249 

1856.5- 3754 
2143.0+0654 



28.06 + 0.07 
26.64 ± 0.03 
27.81 + 0.14 
28.03 ±0.55 
27.17 ±0.11 
25.16 ±0.05 
28.62 ±0.03 



29.81 ±0.06 
28.44 ±0.03 
29.59 ±0.14 
29.81 ±0.55 
28.97 ±0.11 
26.96 ±0.05 
30.37 ±0.03 



29.98 ± 0.06 
28.62 ±0.03 
29.76 ± 0.14 
29.98 ±0.55 
29.15 ± 0.11 
27.13 ±0.05 
30.54 ±0.03 



30.35 ±0.07 
28.96 ± 0.03 
30.11 ±0.14 
30.33 ±0.55 
29.49 ±0.11 
27.48 ± 0.05 
30.91 ±0.03 



30.88 ±0.06 
29.54 ±0.03 
30.67 ±0.14 

30.89 ±0.55 
30.07 ±0.11 
28.05 ±0.05 
31.44 ±0.03 



30.95 ±0.06 
29.60 ±0.03 
30.73 ±0.14 
30.95 ±0.55 
30.13 ±0.11 
28.12 ±0.05 
31.51 ±0.03 



31.66 ± 0.06 
30.32 ±0.03 
31.45 ± 0.14 

31.67 ±0.55 
30.85 ± 0.11 
28.84 ±0.05 
32.22 ± 0.03 



10 B K a 

A cfI a (A) 
Ax/Ay* 



23.16 
2955 
1.93 



25.00 
4520 
1.30 



25.18 
4709 
1.24 



25.50 
5066 
1.58 



26.11 
5844 
0.97 



26.18 
5932 
0.95 



26.91 
7025 
0.79 



Note. — We give the instrument name and filter for all optical observations that we consider; not all objects were 
observed with all combinations, and see Table [3] for the optical observations. The l-<x uncertainties are based on the 
uncertainties in the X-ray spectra (Table [2} and do not include any other systematic terms. 

a The STMAGs for a 10 6 K (86 eV) blackbody with Ay = and normalized to mp45ow = 25. The wavelengths given 
are approximate effective wavelengths for the instrument /filter and a typical INS spectrum. Wavelength-dependent 
extinction A\/Av for a 10 6 K (86 eV) blackbod y , calc ulated for Av = 0.05. Effective wavelengths and extinctions 
follow the definition of van Kcrkwijk & Kulkarni (2001). 



Counterparts and SEDs of Nearby Neutron Stars 



9 



2.5, having the flattest, least Ray leigh- Jeans-like spec- 
trum (as well as the highest excess; see below). No source 
has a slope that is significantly steeper than Rayleigh- 
Jeans. 

Comparing the emission to the X-ray spectrum, there 
is a similarly wide variation in both the amount and slope 
of the optical/UV excess: most sources are not consis- 
tent with the constant excess of RX J1856.5— 3754, and 
the amount of the excess varies significantly, although no 
source has measurements below its X-ray extrapolation. 
At 1500 A the range is from about 3.5 to almost 10, with 
statistically significant variation (x 2 = 134 for 6 DOF) 
even when excluding RX J2143.0+0654 ( X 2 = 46 for 5 
DOF). At longer wavelengths the variation is even larger, 
going from 6 to 50 at 4700 A. Here, RX J2143.0+0654 is 
clearly very different from the rest, but again when ex- 
cluding it there is still significant variation (\ 2 =21 for 
5 DOF). These are beyond deviations possible from the 
X-ray uncertainties, and while X-ray calibration errors, 
differences in X-ray fitting methodology, and variability 
can all contribute to variations in the excess, as noted 
above the shape of the excess is robust. See Figure [T] for 
a full optical-to-X-ray comparison of RX J1856.5— 3754 
and RX J2143.0+0654. 

If the X-ray a nd optical came from different regions 
on the surface ()Braie &: Romanil 120021) , we might ex- 
pect variations in the excess to correlate with changes 
in the X-ray pulsed fraction, as both would be de- 
pendent on geometry: a small X-ray hotspot would 
give rise to both a large optical excess and a large 
pulsed fraction. Some allowances could be made for 
viewing angles, but with 7 sources this would start 
to average out. However, we do not see any such 
correlation. As an example, both RX J1856.5— 3754 
and RX J0420. 0—5022 have similar optical excesses, 
but the p ulsed fraction for RX J 1856.5— 3754 is 
about 1% (|Tiengo fc Mereghettil 120071) while t hat of 
RX J0420.0-5022 is 13% (jHaberl et all l200l . and 
RX J1605. 3+3249 has a larger excess but a small pulsed 
fraction. 

Even so, we can still fit the optical/UV SEDs as part 
of some blackbody, with radius and temperature free 
(Fig. 0}. For RX J1856.5-3754 this is largely degen- 
erate (since it looks like a Rayleigh- Jeans), but for other 
sources the temperature and hence the radius are con- 
strained. Some of these (mostly RX J0420. 0-5022) are 
not likely based on extrapolations to the X-ray band, as 
they would exceed the flux at 100 eV. All fits are gener- 
ally consistent with reduced x 2 ~ lj although there are 
not many degrees of freedom. This fit gives a normal- 
ization of R/d w 400kmkpc _1 for RX J2143.0+0654. 
Since we expect a distance to RX J2143. 0+0654 of 500- 
1000 pc based on the X-ray spectrum, this would imply 
a rather large radius if it is interpreted physically, and 
may only apply in light of the scattering model discussed 
below. Most of the objects fall on a single locus with ra- 
dius increasing as the temperature decreases. Much of 
this just comes from keeping a similar excess among the 
sources (we would expect R oc 1/T), but it is clear that 
RX J1605.3+3249 and RX J2143.0+0654 are inconsis- 
tent with the other objects. 

We saw above that the optical excess does not appear 
to correlate with the X-ray pulsed fraction, suggesting 



100 



\ RX .12143 








i^^X .11,7211 




rx .n;ioS^|^nx .iijsoe 













0.02 0.1 1 

Tow/Tx 

Figure 4. Results of fitting blackbodies to only the optical/UV 
data for the INSs. We show ±1<j (solid) and ±2<r (dashed) limits 
on the temperature and normalization R/d, where in both cases 
we scale the fitted values by the values measured in X-rays (caus- 
ing the distances to cancel). Individual sources are labeled. For 
RX J1856.5-3754 and RX J0420. 0-5022, the results are consis- 
tent with the X-ray temperature at a constant excess, as seen 
from Table [5] Other objects have flatter spectra so the implied 
fits are cooler. For the objects with more than two data points 
and T OU v < T x , RX J0720.4-3125 has \ 2 = 4.3 (2 DOF) and 
RX J1605.3+3249 has x 2 = 5.2 (3 DOF). Extrapolating the op- 
tical fits to higher energies for RX J0420. 0—5022 and to a lesser 
extent RX J1856.5— 3754 might not be consistent with the X-ray 
data. 

that a model like lBraie k. Romanil ([2002) is incomplete. 
Taking the basic spectral parameters {kT along with en- 
ergy of the absorption lines) and the basic rotational pa- 
rameter (pulsed fraction, magnetic field, spin-down lu- 
minosity) we see no clear trends with the optical/UV 
properties that we have measured here. The one possi- 
ble exception is in relating the spectral index to kT: in 
Figure [5] it appears that the hotter objects have smaller 
spectral indices. Much of this correlation comes from 
RX J2143. 0+0654, but even without it the correlation 
appears somew hat significant: S pearman's rank correla- 
tion coefficient (jPress et al.lll992l p. 641) for all objects is 
—0.75 with a null-hypothesis probability of 5%, and this 
drops to —0.6 (21% null hypothesis probability) when 
RX J2143. 0+0654 is excluded. For the small number of 
objects considered here, this is reasonable but not defini- 
tive. 

Such a correlation may imply that the SEDs of all 
of the INSs deviate from a Rayleigh- Jeans tail, but 
that the energy at which they do it depends on their 
temperatures. For the cool RX J1856.5— 3754 and 
RX J0420.0-5022, the optical/UV window still seems 
thermal. But for the hotter RX J1308.6+2127 and 
RX J2143. 0+0654, the thermal portion would be lost 
shortward of 1000 A. The shift does not appear linear, 
as measuring the excess at a wavelength that scales with 
1/kT does not improve the agreement among objects. 

In contrast to some previous attempts to measure the 
optical SED using both ground- and space-based mea- 
surements, we see that all INSs are well fit by a sin- 
gle power-law. Further, our results using only HST 
data are similar to those for RX J1856.5— 3754 and 
RX J0720.4— 3125 that include ground-based data. Our 
conclusions regarding the 5 sources observed here are rel- 
atively insensitive to ^month scale variability, as in many 
(but not all) of the cases the optical and UV observations 
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occurred within weeks of each other, but this is not al- 
ways the case. For instance, RX J1605. 3+3249 now has 
more HST photometry than RX J1856.5— 3754 spanning 
almost 10 years, but still appears smooth. Whether this 
implies that the optical/UV is constant or just that we 
are unlucky is not clear, but future high-quality measure- 
ments should be able to address this. 

4.2. Comparing the INSs to Pulsars and Magnetars 

Prev ious modeling of RX J0720.4-3125 (jKaplan et al.1 
I2003H ) could not distinguish between a single non- 
thermal power-law and a combination of power-laws. 
Unfortunately, with our new data it is still hard to 
make that distinction, as the wide wavelength gap 
between the optical and UV even in the best case 
(RX J1605.3+3249) could easily hide a curved SED. 
Future observations in the 2000-3000 A range might 
help better constrain such spectra and settle this ques- 
tion. We note, though, while the sources with few 
data-points can accommodate a power-law as fiat as 
that for RX J2143. 0+0654 in addition to a steeper, 
thermal component, objects with denser coverage like 
RX J0720.4-3125 and RX J1605. 3+3249 cannot. 

We can set rough upper limits on another power-law 
by requiring that any non-thermal contribution match 
our reddest optical point as well as be below the X- 
ray black body at 1.5 keV, where typically there are no 
more counts (Figure Q]). We also require that such a 
power-law have a slope such that it is falling in pho- 
tons per unit ene rgy toward higher en ergies, so that it 
stays hidden (see lKaplan et al.ll2003bO . T aken together 
with the parallax for RX J1856. 5-3754 ([Walter et al.l 
2010) and assuming all the INSs have similar radii, 
we find non-thermal luminosities integrated over the 2- 
lOkeV band of typically 10 28-30 ergs -1 . At the upper 
end we do not believe that such luminosities could be 
real, as they would be close to 100% of E, in contrast 
to radio pulsars that have non-thermal luminosities of 
~ 10- 3 E in X-rays (jBecker fe Triimperl fl997i ). Even 
without such an extrapolation our optical luminosities 
((0.3 — 1.7) x 10 29 ergs _1 ) are uncomfortably high for a 
typical magnetospheric origin, since it would be > lO -3 ^ 
compared to typical values of ratio s of 10 -7 -10 -6 for ra- 
dio pulsars (jZavlin fc Pa vlov 2004). However, while the 
X-ray emission cannot be powered by E (since X-ray 
luminosities are ;§> E; iKaplan fc van Kerkwiikl [2009a) 
some amount of the optical emission could be related, 
albeit with a different mechanism than that which usu- 
ally operates for pulsars. 

If the deviations from a thermal S ED were from mag- 
netospheric emission as in pulsars (iPavlov et all Il997t 
iShibanov et~aL1 120051: iKargaltsev et al.l l2005f ) ." we might 
expect a correlation in the amount of the excess or the 
optical luminosity with the spin-down luminosity E, as 
this is what drives the X-ray power-laws seen in pulsars 
and may drive the optical too. Again, though, we see no 
such correlation. 

The behavior of the INSs also deviates from what is 
seen in the cases of other relative ly young neutron stars in 
the optical, namely magnetars (Mcrcghctti 20081). Most 
magnetars have only been seen in the infrared (likely 
due to high extinctions). The infrared emission is vari- 
able, and may or may not be linked to the X-ray emis- 
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Figure 5. Optical/UV spectral index (a from Tabled versus X- 
ray temperature kT (from Table [2} for the INSs. The horizontal 
line shows a slope of 4, to be expected for a Rayleigh-Jeans tail. 

sion. The one so urce with confirmed optical emission 
is 4U 014 2+61 (|Hulleman. van Kerkwiik. fc Kulkarnil 
120001 l200l . The optical emission is known to pulse 
at th e X-ray period ()Kern fc Martinll2002t iDhillon et all 
120051). suggesting that it is likely magnetospheric in ori- 
gin (|Zane. Nobili. fc Turollall201lT) . Importantly for this 
compariso n, the optical em ission is not a single power- 
law (jHulleman et al.1 120041 ) . In contrast, the emission 
of the INSs is well fit by a single power-law and is (so 
far) constant in time. The flattening of the SED of 
RX J2143. 0+0654 toward longer wavelengths is some- 
what suggestive of the broadband SED of 4U 0142+61, 
but the other characteristics make it distinct. 

It may be that the similarities with magnetars arise 
from a common origin, which suggests that the optical 
emission results from resonant scattering in the neutron 
stars' magnetospheres [Lvutikov fc Gavriill (|2006l ). How- 
ever, qualitative differences between the INSs and mag- 
netars exist: while the ranges overlap, the magnetars 
typically have higher values of magnetic to thermal en- 
ergy HeB/m e ckT (> 10 4 for the magnetars vs. w 3000 
for the INSs), and the higher fields of the magnetars 
mean that scattering and pair-production happen much 
more quickly. The optical emission of the INSs could 
instead be from an inverted temperature layer (like a 
chromosphere) that is driven by a twisted magnetosphere 
(C. Thompson 2011, pers. comm.). We note that for 
RX J2143. 0+0654, which has the most extreme optical 
excess, we also see that the X-ray pulsations have odd 
harmonics likely indicative of a complex field geometry 
(jKaplan fc van Kerkwiikll2009b[ ). 

4.3. Atmosphere Models for the INSs 

In the context of magnetized atmosphere models for 
neutron stars, the amount of optical/UV excess can de- 
pend on the magnetic field strength. Therefore we ex- 
amined the UV/optical emission predicted by the NS- 
MAX models of NS X-ray spectra constructed for Xspec 
(|Ho et all 120081 ): these models span the range B = 
(1 - 30) x 10 12 G and kT ~ 20 - 400 eV for a partially 
ionized hydrogen atmosphere. We find that, although 
the magnitude of the flux differs from that of a black- 
body, the wavelength dependence still exhibits a A -4 
Rayleigh-Jeans behavior. Deviations from the Rayleigh- 
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Jeans behavior occur as a result of the proton cyclotron 
line at (redshifted) - 2500/B 12 A (with B = 10 12 B 12 G). 
For B < 4 x 10 12 G, the wings of the cyclotron line 
can reproduce the wavelength-behavior seen in Fig. [31 
though not at the magnitude of RX J2143. 0+0654; how- 
ever, it is unclear how strongly this behavior depends 
on temperature (see Fig. [S]). If the deviations from 
Rayleigh- Jeans are due to absorption in the wings of 
the proton cyclotron line, these spectrally-inferred mag- 
netic fie lds are much lower than thos e inferred from 
timing (jKaplan fc van Kerkwiikl I2009al ). A factor of 
two increase in the former may be the result of a 
helium (rather than hydrogen) atmosphere. A possi- 
ble discrepancy between timing and spectral magnetic 
fields may exist when consider i ng the broad X-ray ab- 
sorpt i on features (iHaberll 120071: Ivan Kerkwijk fc Kaplan! 
120071: iKaplan fc van Kerkwiikl l2009b[ ) but there it is in 
the opposite direction (the features are more consistent 
with strong fields), and hence easier to understand in 
terms of higher-order magnetic moments leading to lo- 
cally stronger magnetic fields. We also note that our 
models indicate an increasing flux with increasing wave- 
length similar in magnitude to that seen in the INSs 
(again except for RX J2143.0+0654) for kT ~ 100 eV 
and magnetic fields more similar to what is seen in X- 
ray timing, B w (2 — 3) x 10 13 G; the cause of this is still 
under investigation. 

There exists an uncertainty in the models of strongly- 
magnetized NS atmospheres at long wavelengths. 
The plasma freque ncies for the m odel atmospheres 
constructed in, e.g., IHo et all (|2008l ). lie near the opti- 
cal/UV regime. Classically, emission below the plasma 
frequency should be suppressed, or at least highly 
modified. The suppression would be stronger at longer 
wavelengths (opposite to the behavior seen in Fig. [3]), 
since the plasma frequency increases with density. 
Calculations using an ad-hoc appro ach to account for 
this dense plasma suppression (see IHo et al.l 120031 for 
details) indicate that deviations from Rayleigh- Jeans 
occur at far-UV/soft X-rays (for the likely INS magnetic 
fields B ^> 10 12 G), while at optical wavelengths, the 
flux is suppressed but the wavelength-dependence is not 
affected. Thus the above considerations may still apply. 
Improveme nts on the treatme nt used here are possible 



see e.g., 



2004; 



Brinkmannl 119801: [Turolla. Zane. fc Drake) 
Perez- Azorm. Miralles. & Pons! 120051: 



van Adelsberg et al.l 120051: IHo et al.l 12007) but do 



not qualitatively alter the results, while a fully self- 
consistent model of the emission properties of the 
condensed surface below the atmosphere is beyond the 
scope of this work. 



Regardless of explanation, we have conclusively iden- 
tified optical counterparts to all 7 INSs and shown 
that the relatively straightforward behavior shown by 
RX J1856.5— 3754 is not necessarily the dominant be- 
havior. Future modeling of these sources will have to 
account for this diversity. The quality of the HST data 
has allowed us to measure accurate SEDs and set refer- 
ences for future astrometry. 
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Table 5 

Power-law Fits to HST Photometry 



RX J 




F A (2500 A) 




Excess 




a 


da/dAy 


x 2 /DOF 




(mag) 


(10~ 18 erg cm -2 s _1 A -1 ) 


1500 A 


2500 A 


4700 A 








0420.0-5022 


0.11 


0.34 + 0.03 


9.6 + 1.3 


8.3 + 0.8 


7.2 + 0.7 


4.20 + 0.22 


1.2 


0/0 


0720.0-3125 


0.06 


0.92 ± 0.05 


4.6 + 0.3 


6.1 + 0.4 


8.6 + 0.5 


3.43 + 0.12 


1.4 


1.5/2 


0806.4-4123 


0.10 


0.22 + 0.02 


3.6 + 0.6 


4.3 + 0.7 


5.5 + 0.9 


3.62 + 0.20 


1.2 


0/0 


1308.6+2127 


0.10 


0.19 + 0.01 


3.8 + 2.1 


4.5 ± 2.5 


5.6 + 3.2 


3.63 + 0.14 


1.2 


0.4/1 


1605.3+3249 


0.05 


0.68 ± 0.02 


5.0 + 0.5 


7.2 + 0.8 


11.6 + 1.3 


3.23 + 0.07 


1.3 


3.5/3 


1856.5-3754 


0.05 


4.25 + 0.07 


7.1 + 0.4 


7.1 + 0.3 


7.3 + 0.3 


3.94 ± 0.06 


1.9 


2.6/2 


2143.0+0654 


0.13 


0.48 ± 0.02 


9.5 + 0.4 


19.9 + 0.9 


49.8 + 2.3 


2.53 + 0.08 


1.2 


0/0 



Note. — Wc give the nominal extinction Ay^i (computed from A^h in Tabic [2j, the absorbed flux density at 2500 A, the 
power-law index a (defined such that i*\ oc A~ a ). and the \ 2 and degrees of freedom for the fit. We also list the excesses 
over the X-ray blackbody at 2500 A (where the normalization is calculated) and at 1500 A (appropriate for F140LP) and 4700 A 
(appropriate for F475W). Uncertainties on the X-ray fits are included in the given uncertainties on the excess. Extinction at 
2500 A is taken as j ^ 2 sooa^ v = ^c a l so gi ye the change in power-law index with changes in extinction, da/dAy, i.e., 

a(Ay) = a(Ay Ic i) + da/dAy(Ay - A VlX et)- 
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